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Light emission in sonoluminescence is shown to be a lasing process with a wide gain bandwidth. 
Population inversion of the gas molecules inside the bubble is achieved by hydrodynamical pumping. 
Analytic expressions are derived for the sonoluminescence pulse: shape, pulsewidth, and their code- 
pendence on the spectrum and intensity in physically relevant regimes. A detailed comparison with 
experiments(R. Hiller et ai, Phys. Rev. Lett. 80, 1090 (1998); M.J. Moran et ai, Phys. Rev. Lett. 
80, 4987 (1998); R. Pecha et ai, Phys. Rev. Lett. 81, 717 (1998)) suggests excellent agreement. 
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The phenomenon of sonoluminescence (SL) emission 
of a short pulse of light in an acoustically levitated bubble 
during its final stage of compression — has attracted a lot 
of interest jl| . Unfortunately, there does not exist a com- 
plete theory that formally explains all aspects of the emit- 
ted light §-|. Light emission in SL can be understood 
without invoking bubble dynamics in contrast to the re- 
cent claims |j) that a modified bubble dynamics theory 
could numerically account for the nature of the pulse 
|^]-|5| . We argue that population inversion in the medium 
of gas molecules, which must contain metastable states, 
generates the SL pulse; towards that end, we explicitly 
calculate measurable pulse characteristics and compare 
with the experimental observations However, bub- 

ble dynamics is important in evaluating the temperature 
inside the bubble with precision from ambient parame- 
ters such as pressure, bath temperature and equilibrium 
radius. 

We argue that light emission in SL results |7| mainly 
due to the following two steps: (i) excitation of both gas 
atoms and molecules containing metastable states at a 
rate high enough to sustain population inversion against 
dephasing and loss, (ii) spontaneous and stimulated emis- 
sion of the excited atomic/molecular levels, collective due 
to a high density. We consider a realistic model of an ar- 
gon gas bubble and evaluate various parameters of the 
emitted light pulse in view of the proposed mechanism 
of sonoluminescence as a broadband laser |Q — though an 
inefficient one due to the presence of strong dephasing 
and loss. In this exercise of SL from an argon bubble @, 
we derive explicit analytic expressions for (i) the tem- 
poral shape of the pulse in various relevant regimes, (ii) 
the pulsewidth consisting of the rise time — denoting the 
growth of the pulse, and the fall time — denoting the de- 
cay, both compared with recent streak camera experi- 
ments in ref. |@J|, (iii) their individual dependence on 
wavelength, compared with experiments in ref. , (iv) 
the intensity of the emitted light and its dependence on 
the pulsewidth, compared with the experiments in ref. Q , 
and (v) the spectrum at different ranges of wavelength, 
compared with experiments in ref. . Our results are 
in excellent agreement with these experiments 

Chemical Requirement: It is known that the gas tem- 
perature inside an argon gas bubble can be as high as 



10,000 K to 40,000 K (0.86-3.5 eV) @ due to hydro- 
dynamic collapse 0. Heating to such a high tempera- 
ture is sufficient to cause electronic excitation of argon 
gas atoms. This process is necessary in creating popula- 
tion inversion in atoms with metastable states essential to 
standard lasing process. However, for hotter bubbles of 
argon and other rare gases fll|l with higher atomic num- 
ber (smaller ionization potential) electronically bound 
excited molecular states are formed, though the inter- 
action between the constituent atoms is mostly repulsive 
Jl2| | : Ar + Ar* — * Ar*. , which then decay to ground state 
by emitting a photon: Ar*, — » Ar+Ar+hu. These excimer 
states are dissociative in the electronic ground state. 
Spontaneous and stimulated decay of such molecular ex- 
cited states describe the underlying physical mechanism 
of excimer laser systems such as Ar*, or Xe?, lasers. Ex- 
cimers with strongly repulsive ground states emit broad 
continuum radiation whereas the weakly bound excimers 
typically display a radiation spectrum with rotational 
and vibrational structures, which gets homogeneously 
broadened at high pressure. The binding energy of the 
excimers depends on potential-energy curves which are 
extremely complex, including families of closely nested 
curves of excited states as well as level crossings linking 
specific levels. After excitation, the system evolves from 
the top to the bottom of these curves by relaxing through 
a nested sequence of electronic states. The chemical ki- 
netic ordering originating from the particle interactions 
of rare gases is such that no direct thermal path is avail- 
able, and the only accessible decay path is radiation. We 
argue that this kind of radiation is what gives the light 
in sonoluminescence. 

Let us now consider the chemical kinetics of argon in- 
side the bubble specifically for this excercise. During the 
compression of the bubble at high pressures (P > 10 3 
atm) and temperatures (ksT e ~ 3.5 eV), argon is excited 
(e~~ + Ar — > Ar* + e~) or ionized (e~ + Ar — > Ar + + 2e~) 
with a population of e~ Eian / kBT or e~ Eex ^ kBT reS p ec _ 
tively. Since the ionization energy of argon is ~ 12 
eV, this corresponds to a small degree (~ 3%) of ion- 
ization at a temperature in the range of 0.86-3.5 eV. 
The atomic excitation to Ar* is much more efficient. As 
we will see, a large amount of initial ionization is not 
important for the formation of excimer states. Three- 
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body recombination reactions rapidly dimerize the 
excited atoms (Ar* + 2Ar — > Ar*, + Ar) or ionized states 
(Ar+ + 2Ar -> ArJ + Ar) @. Ar gon excimers formed 
mostly in many excited states quench to lower states at 
high pressures p3] and the excitation resides mostly in 
the lowest excited levels of atomic and molecular levels. 
The radiative lifetime T\ of the most relevant excited 
state Ar?, (for the l u molecular state which dominates 
over the 0+ state) is known Q to be 3 x 10~ 6 s. 

Presence of a small amount (<C 10%) of a heavier rare 
gas in a large amount of a lighter rare gas || results in 
enhanced ionization. These impurities can dimerize, re- 
combine and react to form various diatomic gas excimers. 
If the excited species are sufficiently energetic, then they 
react with the heavier rare gas atoms, Xe in Ar for exam- 
ple, to form ions in Penning (Ar* + Xe — > Ar + Xe + + e~ ) 
or associative (Ar*+Xe — > ArXe + + e~) ionization; these 
reactions are very fast and occur at every collision. If 
Penning or associative ionization is not energetically fa- 
vored, then the energy of the lighter species transfers to 
the heavier species by reactions Ar* + Xe — > Ar + Xe* , 
and Ar* + Xe -> 2Ar + Xe* 

This roughly describes the processes through which ar- 
gon gas excimers are formed inside a sonoluminescing 
bubble. In the same spirit, one anticipates other rare gas, 
rare-gas-halogen excimers and quasi-bound metal vapor 
excimers Jl2| to show lasing in the SL configuration. The 
above chemical kinetics analysis applies to all other kinds 
of excimers, though the details of the collision mecha- 
nisms and reaction rates are different. In the following 
we discuss the second step of the process which gives rise 
to the light emission in SL from an argon bubble: decay 
of the excimer states Ar*,. 

Hamiltonian: Let us define a transition operator Rf b 
from any state |o) to \b) of the i-th molecule: Rf b \a)i — 
\b)i. The hamiltonian of the molecular system compris- 
ing of multiple levels can then be described by H at = 
J^i^H^o + u e )R1 e + J^i.gh^gRi 9 ') where the index % 
runs over all the molecules, and the index e and g de- 
fine the energy of a given excited state |e)j and ground 
state \g)i with respect to an average energy %u)q. Tiloq 
is chosen, without any loss of generality, to be the dif- 
ference in energy between the lowest excited state and 
the lowest ground state. The hamiltonian of the field 
is given by the sum over all modes k and polarization 
a such that H ra d = J2ka ^k^i^k- The interaction 
term between t he mo lecular states and the field modes is 
Hint = Eke 34 \J l^v[ a KRTd ge {6 a ■ e ge ) + h.c.}. 

Polarization and Population Inversion: Collective 
dipole moment of the molecules can be expressed in terms 
of a polarization operator, constructed from the molecu- 
lar raising and lowering operators Rf(u>). Denoting the 
projection of the corresponding molecular dipole matrix 
element d ge along the polarization direction as d u , we 
define the following: P ± (r,uj,t) = J2i d u Rf(r — r i; u>, t). 

N(r,to,t) = J2i R-i 3 \ r ~ Ti^'t) is the number of excited 
molecules at a transition frequency w. 



Maxwell- Block Equations: The common electric field 
created by the excited molecules is related to the polar- 
ization by the Maxwell's equation: 



d 2 E+(r,t) 
dt 2 



c 2 V x VE+(r,<) = - 



1 <9 2 p-(r,i) 



dt 2 



(1) 



From the Heisenberg's equation of motion for N and P: 



dN _ i 
~dt ~ti 



(5P+ n] 2 

— = ilu P + + 2i—e a ■ %{E-N + NE+)], (2) 

which define the Maxwell-Bloch equations (MBE) in op- 
erator form [p^[ . Electric field and polarization in an 
arbitrary direction, the direction of the detector z for ex- 
ample, can be expressed in the slow- varying approxima- 
tion in a complex form as E + (r, f) = z£ + (z,t)e l ^ ot ~ kr ^ 
and P+(r,w,t) = zV + {z, w, t)e^ ot ^ kr \ where t is re- 
placed by the retarded time t — r/c, introduced for con- 
venience, and £ and V are the slow- varying envelopes. 
To describe the experiments it is essential to include the 
effect of dephasing and loss p9[ : k denotes propagation 
loss of the electric field and is related to the refraction 
index rj by kc/u; = r\ — in, I/T2 is the dephasing rate of 
polarization, 1/T\ is the spontaneous emission rate, Ap 
is a noise source of polarization due to the background 
field, and finally, A describes the population inversion by 
the excitation due to hydrodynamic pumping. Eqs. (1) 
and (2) now take the following form: 



d£+ 

dz 

cVP+(r,uj,t) 



= — k£ 



2e c 



V (r,t,u)g(uj)du. 



A p - (i- - ikv)V + + ^-£-M{r,w,t), 

= A - fe^ - l -(v+£+ £-V-). (3) 
dt Ti h 

The first equation describes the electric field created by 
the polarization of the molecular excitations of the rele- 
vant manifold of states minus the loss of the field due to 
propagation. A p represents the noise polarization due to 
population inversion and the blackbody radiation, which 
act as an equivalent input field to stimulate the decay. 
The effect of Doppler broadening relevant at high tem- 
peratures is represented by the term ikvV + . In fact, it is 
not necessary to account for all these phenomenological 
factors to describe the essential features of the recent SL 
experiments. To that end, we shall obtain the equation 
for the radiated pulse for a simplified case in the time do- 
main, by replacing the opeartors £, V and Af in the MBE 
by their complex functional forms. From the first MBE 
and its complex conjugate, after multiplying by ccqA£ *, 
an expression for the radiated intensity I is found: 

81/ dz = -2kI + {<j A/2)2Im[£ [ V(r,w,t)du], (4) 
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where A is the area of the cross-section perpendicular 
to the direction z. J Vdw gives a polarization at a cen- 
tral frequency luq averaged over an effective bandwidth 
(Auj) e ff. Differentiating with respect to t and using 
Eq. (3) one obtains: 



81 

dt 



= -(1 + 2k£) 



T-2 



2uj Ad 2 NL 

h 



I + uoALApE, (5) 



where L is a characteristic length scale of the active 
medium. The first term on the right hand side de- 
scribes dephasing or inhomogeneous broadening (I/T2) 
and propagation loss (k). The second term denotes the 
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I. The third term 



source, and it can be written as -^r T 
is the polarization source for the input field that drives 
the decay: 
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" Lo AL TR T* ff 



n 2 



n e Huj /k B T _ 2 J 



(6) 



The effective dephasing time for a gain medium 



is t: 



eff 



T 2 (2a L/7r) 1 / 2 , where a Q L 



T 2 /T R is 

-1 



the gain over the length L, and T r (u)q, Au> e f f 
N(ivo,Auj e ff)/Ti is the collective decay rate at ujq. The 
second term gives polarization due to blackbody radia- 
tion interacting within Aw e // = l/Tj , and the first 
term inside the bracket is the excitation spectrum of the 
population inversion g{u]Q, Ao; e //) = n 2 /(n 2 — In 
thermal equilibrium n 2 — nie ~ huJ o/k B T ^ ^hus the right 
hand side of Eq. (6) vanishes, implying the absence of 
population inversion. 

Coupled equations for intensity and population are ob- 
tained from Eq. (5) and (6): 



-(1 + 2kL] 



dt 

^jj[g(ujo,Auj eff ) 



TrT 2 
dN 
~dt 



N 



AAL 



I 3A^c NL 
1 ^ j 

T 2 4tt 2 Ti 

I 1 

e huj /k B T _ ]_J' 
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hujQ 



(1 + 2kL), 



(7) 



(8) 



T[ = T 2 if T 2 < T x as in the case of SL, else T[ = T x . 

Pulse Shape in an Ideal Case: In an idealistic case 
one neglects dephasing, T 1 ~ 1 ,T 2 ~ 1 = 0, k = 0, and 
A, A p = 0, retaining the effects in the initial condi- 
tions. We define a time- and space-dependent tipping 
angle 9(z, t) in terms of the population and polariza- 
tion in the medium so that TV = uq cos#(z, t), and V = 
idn sin 9 (z, t) e 1 ^ . The Bloch equations at a given u> re- 
duce to (2d/K)£ — e^dO/dt, and the Maxwell's equation 
to d 2 9/dzdt = (TN/L^smd, where FN = T X N = T^ 1 , 
and Lff, is the length scale over which t he p olarization 
contributes to the electric field coherently pQ| . The MBE 
can be written as 

9"(u) + 6'{u)/u-sm9 = 0, (9) 

in a dimensionless variable u — 2(^-t£-) 1 ^ 2 , demonstrat- 
ing clearly the natural emergence of a collective timescale 



T R . Solutions to the damped pendulum equations of this 
type, known as sine-gordon equations, include a wide 
range of elliptic functions with the constraint 9(0) = 9 , 
the input noise field, and 6'(0) = 0. If L -C L^, then the 
entire active medium is coherent and one obtains pure 
superradiance. Solution of the equation in that limit 
with the first term omitted yields an expression for in- 
tensity / cx sech 2 (t - T D /T R ), where T D = T R ln(N) 
is the pulse rise time Q. When z > L^, then one ex- 
pects the occurence of several maxima(ringing) in inten- 
sity. The field emitted by molecules in the central co- 
herent volume(z < L^) are reabsorbed by the molecules 
outside this regime (L^ < z < R m in) to give stimulated 
emission. Hints of this behavior may have been experi- 
mentally observed jf|. 

Pulse Rise Time: The pulse rise time is defined as the 
time at which the first peak is observed from the time 
when the sample is completely inverted. In the small 
angle limit (sinO ~ 9), Eq. ^ has the standard Bessel 
function solution of the linearized sine-gordon equation: 
9(q) — 9(0)Jo(iq). Using the asymptotic limit, Jo(iq) ~ 
e q / ^/2/Kq, and requiring 9(q,t = Td) — 1 for the maxi- 
mum reached at time Td, one obtains a transcedental ex- 
pression: T D = (T R /4)[ln9{0) - (l/4)ln(T R /16TT 2 T D )} 2 ■ 
The second term inside the bracket is negligible com- 
pared to the first for the experimentally relevant pa- 
rameters, and T R /l6ir 2 T D < 10. Using Eq. M) one 
finds 9{0) = [2d/h)T R £ mput cx (g(u> ) + ( e ^^ T - 
1) . Thus in the range where the excitation spec- 
trum is thermal, ((e^AsT _ ^-1 > 5 ( Wq )) ; t d = 
(T R /4)[ln(k B T/hu))} 2 - constant, for k B T > hoj , and 
T D = (T R /4)(huj/2k B T) 2 - constant, for k B T < huj - In 
the nonthermal regime of the excitation spectrum, the 
rise time is independent of temperature. In experiments, 
since T is usually much larger than hu>/k B (~ 2500-11000 
K in the visible range) , dependence of Td on A and T is 
expected to be logarithmically slow, or almost absent, in 
agreement with the streak camera measurements . 

Pulse Fall Time: In the limit where dephasing by in- 
homogeneous broadening becomes dominant, T 2 <C T R , 
an expression for the radiated intensity can be obtained 
by retaining all the terms in Eq. (|?|): 



I(t) = /(0)e- { / Ti + n 



hujQ 



+ 



TrT 2 
1 



e huj /k B T _ 1 



eff L5(^o,Aw e// ) 

}(l-e- t/T >), (10) 



where the effective fall time is given by r i 1 = T 2 1 (1 + 
2kL) — 3\qcL(NT^ 1 ) / 4n 2 . In the absence of propagation 
loss (k ~ 0)), and negligible superradiance rate 1/Tr, the 
fall time of the pulse is dominated by dephasing due to 
inhomogeneous broadening, Ti — T 2 . An expression for 
N is similarly found to be 

N(t) = N(0)e~ t/T2 + T[KAL{1 - e^ 2 ) 

_M?)r 1 '(i + 2 K L)(i-r t / T! ) (11) 
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as T[ ~ T 2 « Ti. Thus in the dephasing dominated 
regime both I(t) and N(t) are seen to decay exponentially 
with a characteristic time of T 2 rather than with a sech 
form which is expected in the regime where dephasing is 
negligible, T 2 3> Tr. The asymmetry in the pulse shape 
is naturally anticipated, since To 7^ Tg. The presence 
of the two timescales and their somewhat independent 
evolution is clearly manifest in the asymmetry observed 
in streak camera experiments(See Figs. 3 and 5 of ref. ||, 
Fig. 4 of ref. Q.) along with the expected exponential 
decay of the pulse. 

What determines the dephasing time T2? Normally at 
high temperatures T 2 is determined by Doppler broaden- 
ing Jig . The fractional Doppler broadening, Auid/w = 
^{8ln2)k B T/Mc 2 , where Mc 2 is the rest mass energy of 
the molecule. Defining the inhomogeneous lifetime 
as T 2 = 3/Aw rf , one obtains, for Ar^ at T= 20000 K, 
and A = 400 nm, Afd = 3.75 GHz, and a correspond- 
ing T2 = 130 ps, comparable to the typical measured fall 
time in the experiments [0-H]. Note that T 2 oc 1/VT and 
T2 cx A, denoting that Doppler dephasing time is larger 
for longer wavelengths, which may explain why the rise 
time is observed to be almost independent of A whereas 
the fall time increases linearly with A as observed in Fig. 3 
of ref. [B. But at high pressures, there is a transition 
from Doppler broadening to collisional processes, which 
depend on the pressure as well as the collision partner. 
Correlation experiments described in ref. || may be in 
this regime. 

Energy emitted coherently can be obtained by inte- 
grating Eq.(|o]): /(*) = Nhu/T R sin 2 0e~ t/T '\ the maxi- 
mum radiated energy is found to be NTiuj(T 2 /Tji), indi- 
cating that only a fraction of the energy is emitted coher- 
ently, only for a time T2 shorter than Tr. The rest of the 
energy is emitted incoherently Q (T2/TR = (aoL) <C 
1). Thus in the dephasing limited regime one finds that 
the emitted intensity at first linearly increases with in- 
creasing pulsewidth, saturating at a value Tr beyond 
which T2 becomes irrelevant. The radiated pulsewidth 
is determined only by Tr, as observed in the experi- 
ments in ref. ||(See for instance Fig. 2 ). Tr can be 
estimated as min{TR(o;o)} — T\/N{ijJq). A typical esti- 
mateis Tx/Nfao, Auj eff ) = (Ti/N)(f max - f mm )T 2 . The 
total N of 72 x 10 7 is obtained at a density of 100 times 
that of STP for a volume of (tt/6)(A = 800nm) 3 . For the 
visible range, this gives a Tr of 600 ps for T\ = 3 x 10 _6 s 



for Ar 2 (L4J and an average T 2 of 130 ps, consistent with 
experiments ||. 

Emitted Spectrum: The energy emitted at a given 
u>o can be obtained from Eq. (|l0|). Total energy 
emitted is found by integrating the energy density 
(ira Q L/2) l l 2 \g(oj Q ) + ( e W fc fl T _ OV e r the fre- 

quency range with the appropriate density of states. At 
longer wavelengths (or eqivalently at higher tempera- 
tures) ( e huJ o/k B T _ 1 3> g(u)o), thus the emission spec- 
trum is almost thermal or blackbody-like. This explains 
why at longer wavelengths (or equivalcntly in hotter bub- 



bles) all the SL spectrums are similar P,||,p[ . This is not 
a surprise since the emitted light in this case is actually 
an amplification (with a gain of a^L) of the input black- 
body stimulation. The molecular electronic structure is 
evident at shorter wavelengths when the blackbody term 
is negligible compared to g(cJo)- 

Dependence of the pulse characteristics on the bath 
temperature and ambient pressure can be studied in our 
model by understanding how they affect the gas temper- 
ature T, the intrinsic parameter in our theory, for which 
the reader is referred to the work in ref. . 

In conclusion, our model of sonoluminescence as a las- 
ing process with a broadband optical gain explains the 
experimentally observed pulse characteristics, the spec- 
trum as well as the parametric dependences of the emit- 
ted light analytically. Our calculations and comparisons 
with experiments will be presented in greater detail else- 
where. I thank S.V. Khare and R. Lifshitz for many 
valuable comments. I acknowledge the formal support of 
Michael Roukes. 
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